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Abstract 

Advances in the amplification and manipulation of ultrashort laser pulses has led 
to revolutions in several areas. Examples include chirped pulse amplification for gen¬ 
erating high peak-power lasers, power-scalable amplification techniques, pulse shaping 
via modulation of spatially-dispersed laser pulses, and efficient frequency-mixing in 
quasi-phase-matched nonlinear crystals to access new spectral regions. In this work, we 
introduce and demonstrate a new platform for nonlinear optics which has the potential 
to combine all of these separate functionalities (pulse amplification, frequency transfer, 
and pulse shaping) into a single monolithic device. Moreover, our approach simulta¬ 
neously offers solutions to the performance-limiting issues in the conventionally-used 


1 



techniques, and supports scaling in power and bandwidth of the laser source. The ap¬ 
proach is based on two-dimensional patterning of quasi-phase-matching gratings com¬ 
bined with optical parametric interactions involving spatially dispersed laser pulses. 
Our proof of principle experiment demonstrates this new paradigm via mid-infrared 
optical parametric chirped pulse amplification of few-cycle pulses. 


1 Introduction 

Intense nltrashort laser pulses play a pivotal role in numerous areas of science and tech¬ 
nology. In industry they have enabled advances in micromachining, while in science they 
enable a broad class of intense light-matter interactions, with applications such as resolving 
attosecond dynamics in atoms and molecules , generation of soft x-ray radiation via high 
harmonic generation [^|^, and driving relativistic laser-plasma processes |^. There is thus 
major interest in advancing these laser sources on several fronts, including wavelength (to¬ 
wards the mid-infrared), pulse duration (towards single-optical-cycle pulses), and repetition 
rate (to rapidly perform experiments and obtain good signal to noise ratios, while avoiding 
detrimental high-intensity effects). A compelling approach to generate the required sources 
is optical parametric chirped pulse amplihcation (OPCPA) [^, where an intense and narrow- 
bandwidth pump pulse amplihes a broad-band, temporally chirped signal pulse in a nonlinear 
crystal. In recent years there has been rapid growth in high-power ultrafast lasers (7-11 


With OPCPA, the energy of these lasers can be transferred to few-optical-cycle pulses with 
user-chosen center wavelengths from the visible to far-infrared. 

Nonetheless, optical parametric processes, including but not limited to OPCPA, present 
numerous challenges. Complicated non-collinear beam geometries and high laser intensities 
close to the damage threshold are needed to achieve phase-matched amplihcation over an 
ultra-broadband spectrum. Moreover, amplihcation occurs at all points in space and time 
of the pump pulse, and exhibits back-conversion of energy from the signal to the pump if 


the intensity is too high 12 . Thus, maintaining the desired interaction across the spatial 
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and temporal/spectral profiles of the interacting ultrashort waveforms remains a demanding 
problem. 

Recently, a complementary approach to OP CPA has been introduced, termed frequency 
domain optical parametric amplihcation (FOPA) [T3p4]. The seed-pulse is dispersed spatially 


via a 4-f arrangement analogous to a pulse shaper 15 , with amplihcation occurring at the 
Fourier plane. By hlling this plane with several birefringent phase-matching crystals placed 


side by side, as in the hrst experimental demonstration of the FOPA technique 14 , the 
phase-matching condition for different spectral regions can be adjusted separately, thereby 
relaxing one of the key constraints of conventional OPCPA systems. Moreover, because the 
seed is spatially chirped, its effective pulse duration can be matched to that of the few-ps 
pump pulse, allowing for efficient energy transfer. On the other hand, drawbacks to this 
powerful approach are that the optical path length through each crystal must be precisely 
matched, the complexity scales with the number of crystals used, and pre-pulses can be 
introduced by any diffraction from the edges of the crystals. 

Here, we introduce and demonstrate a new paradigm for nonlinear-optical devices based 
on combining spatially dispersed laser pulses with a two-dimensionally patterned quasi-phase- 
matching (QPM) medium. This approach represents an extremely versatile yet experimen¬ 
tally simple platform, allowing for the limitations of existing nonlinear devices, including 
the above-mentioned drawbacks of the FOPA, to be overcome systematically by lithographic 
patterning of optimal 2D-QPM gratings. We experimentally demonstrate the approach with 
a mid-infrared FOPA. 


2 Two-dimensional quasi-phase-matching devices 


In QPM 12,17, the sign of the nonlinear coefficient is periodically or aperiodically inverted. 


augmenting the phase-matching condition with a term Kg to yield |kp — kg — ki — Kg| 0, 
where kj are the wavevectors of the interacting waves. In periodically poled ferroelectric 
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Figure 1: Illustration of 2D-QPM concept for frequency-domain optical parametric pro¬ 
cesses. (a) Schematic of the relevant experimental conhgurations, with a spatially chirped 
ultrashort pulse incident on a two dimensionally patterned QPM (2D-QPM) medium de¬ 
signed to individually address the different spectral components of the pulse, (b) Signal 
wavelength versus transverse position in a 4-f pulse shaper (grating frequency 75 lines/mm, 
f=200 mm), and QPM period assuming a 1064-nm pump [^. (c) Example pump intensity 
prohle across transverse position, and corresponding effective length to achieve a flat small- 
signal gain prohle. (d) QPM period (in /im) obtained by combining (b) and (c), together 
with a smooth variation of the period along the longitudinal direction (see Methods), (e) A 
selection of ferroelectric domain prohles (every 1/20*^ domain) for the QPM period mapping 
described by (d). (f) Absolute QPM phase 4>qpm at the input position z = 0 correspond¬ 
ing to part (e). (g) Contours of constant period corresponding to part (d). (h) Stretched 
image of the fabricated QPM grating in the 1-mm-thick MgOiLiNbOs crystal used for our 
experiments. The image was constructed via a series of microscope images of the +z facet of 
the crystal along the transverse direction. The surface was etched to reveal the ferroelectric 
domain inversions. 


materials such as LiNbOa, a lithography mask dehnes the QPM grating with high robust¬ 


ness 


17 19 . Thus, whereas birefringent phase-matching relies only on favorable material 


properties, QPM media can be freely engineered via lithography. For example, chirped QPM 
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gratings can extend the phase-matching bandwidth well beyond that of periodic QPM grat¬ 
ings [20|^ . Here we take a much more general approach by using fully two-dimensional QPM 
(2D-QPM) patterns to tailor the nonlinear interactions experienced by spatially-separated 
spectral components, as illustrated conceptually in Figure [^a). 

First, the QPM period can be varied in the transverse direction such that each spectral 
component is perfectly phase-matched. Figure [^b) shows an example of this procedure for 
a mid-infrared pulse. In general, we can vary the QPM period continuously according to 
the exact trajectory required by material dispersion, with no inherent bandwidth constraint. 
Moreover, the linear properties of the QPM crystal remain homogeneous, so only a single, 
monolithic, plane-parallel crystal is required. 

Next, hgure [^c) shows an example Gaussian intensity prohle of the pump beam. Nor¬ 
mally, such changes in pump intensity would signihcantly change the gain for different spectral 
components. However, the 2D-QPM concept enables variation of the QPM grating proper¬ 
ties along the beam propagation direction as well. As an example to illustrate this general 
capability, we show in hgure [^c) how the effective length L can be matched to the pump 
beam’s intensity prohle, thereby modifying the nonlinear interaction in order to hatten the 
small-signal gain spectrum. 

To design a practical QPM grating prohle meeting this criterion, we introduce a smooth, 
nonlinear variation in the QPM period. A corresponding map of QPM period is shown in 
hgure [^d). To fabricate this design, we introduce the absolute phase of the QPM structure, 
(Pqpm{x,z), given by 

(j)QPM{x,z) = (j)QPM{x,0)+ Kg{x,z')dz', (1) 

Jo 

where the local QPM period is 27i/Kg, and the input QPM phase (j)QPM{x,0) is a design 
degree of freedom (see Methods). Given (j)QPM, the nonlinear coefficient satishes d{x,z) = 
sgn.{cos{(j)QPM{x, z))). The QPM phase is imparted to any waves generated during the non¬ 
linear process: therefore, it introduces an opportunity to create an arbitrary phase mask for 
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pulse shaping purposes. We discuss this potential in section 

The chosen (j)QPM{x, 0) function is shown in hgurej^f), while hgure[^e) shows several of 
the resulting ferroelectric domains. We emphasize that the domains have no discontinuities, 
but have a signihcant curvature. This unique feature strongly contrasts with conventional 
mechanically-tunable QPM devices, which utilize multiple separate gratings for discrete tun¬ 
ing or straight but “fanned” ferroelectric domains in fan-out devices used for continuous 
tuning ( 2 ^ . 

To conhrm that these domain prohles could be fabricated, we inspected the entire width 
of the devices used, as illustrated in hgure [^h), where we show a part of the successfully 
fabricated grating. The hgure is stretched along the longitudinal direction to make the 
individual curved domains visible. There were no noticeable domain errors in the entire 
1-mm-thick MgOTiNbOs crystal (dimensions: 25-mm width by 12-mm length). 

3 Frequency-domain parametric interactions in 2D-QPM 
media 

We next elucidate the parametric process occurring in 2D-QPM media. We focus on 4-f pulse 
shaper arrangements as the means to introduce spatial chirp. Figure [^a) shows a simulation 
of a simplihed situation involving a continuous-wave pump and signal, but including the 
longitudinal variation of the grating. Exponential amplihcation occurs, followed by depletion 
of the pump, followed by a rapid change of the QPM period to “turn off” the interaction. 

The capability to turn off or modify the parametric interaction in other ways within the 
nonlinear crystal is unique to structured QPM devices. Moreover, we show in hgure |^b) how 
these modihcations can be performed in a frequency-dependent fashion. The hgure shows 
a series of the simulations from hgure |^a) for diherent transverse positions: by turning oh 
the interaction after the relevant ehective length, back-conversion of energy to the pump is 
suppressed. Importantly, although frequency-dependent modihcations can be accomplished 
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Figure 2: Modeling of frequency domain optical parametric amplification (FOPA) in a 2D- 
QPM medium, (a) Plane- and continuous-wave interaction in a longitudinally-varying QPM 
grating, showing the procedure used to switch off the parametric amplihcation after a certain 
distance through the crystal, (b) A series of simulations like (a), showing the evolution of 
the pump as a function of transverse position. At each transverse position, we perform a 
separate plane- and continuous wave simulation, with the pump intensity and effective length 
according to hgurej^c). (c) Full spatiotemporal simulation of the FOPA process. The hgure 
shows the output electric field envelopes of the pump and signal for the transverse position 
X = 0. (d) Evolution of the pump fluence through the crystal as a function of transverse 
position, (e) Output signal spectra for three cases: the 2D-QPM grating pattern introduced 
here; a simpler fanout pattern with no longitudinal variation; and the simplest case of a 
standard periodic grating, (f) Group delay spectra for the signal and idler, assuming the 
2D-QPM grating pattern. For case 2, the pattern is flipped with respect to the longitudinal 
coordinate, changing the phase mask seen by the idler wave (derived in Supplementary section 
3) but not changing the gain. 


in longitudinally-chirped QPM devices using one or multiple QPM gratings 20,24,25 , such 


devices are ultimately limited by coupling between different parts of the spectrum 26 . In 
contrast, the spatial chirp in the FOPA decouples the spectral components more robustly, 
enabling greater flexibility. 

To reveal the extent of this decoupling for the more subtle case involving real pulsed 
beams, the complete spatiotemporal profile of the signal must be considered. Our analysis 
in Supplementary section 2, where we derive this profile, reveals a correspondence between 
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the spatial profile at the input diffraction grating of the 4-f setup and the temporal prohle in 
the Fourier plane. Consequently, the duration of the spatially-chirped signal pulse is 




(A) 


A Ax Win{X) 
cAA / 


( 2 ) 


where AA is the range of wavelengths involved, and Ax is the spatial extent of those wave¬ 
lengths. As well as highlighting an important physical aspect of frequency-domain nonlinear 
optics, equation ([^ allows comparison of Tg// to the pump duration, which yields design 
guidelines for efficient operation. 

To capture the nonlinear dynamics of the amplihcation process, we developed a numerical 
model for nonlinear mixing processes in 2D-QPM media (see Methods). Figurej^c) shows the 
input and output of a full spatiotemporal simulation, accounting for propagation coordinate 
2 ;, transverse dimension x, and time t (2-l-lD). The input pump (dashed red) has duration 
14 ps (FWHM). The effective input signal (dashed blue) has duration 5.2 ps (FWHM). 
The regions of the pump overlapped with the signal experience strong depletion, but since 
the signal is shorter than the pump in this example, the temporal wings of the pump remain 
undepleted. Figure |^d) shows how this issue manifests as a function of transverse beam 
position, by integrating over the time coordinate. In contrast to|^b), complete depletion of 
the pump does not occur. 

The importance of a fully two-dimensional pattern is shown in hgure |^e), which plots 
the simulated output spectrum for three cases. The solid blue curve corresponds to a 2D- 
QPM pattern [see hgure [^d)], showing amplihcation of the full input spectrum. The red 
curve uses a “fanout” grating [QPM period varied transverse to the beam according to hgure 
l^b)] but has no longitudinal variation. There is a substantial reduction of the bandwidth, 
and over-driving the device to recover the lost bandwidth introduces strong spatiotemporal 
distortions. Finally, the black curve models a simple periodic grating: in this case, there is a 
drastic reduction in bandwidth. 

Beyond amplifying the interacting waves, the QPM device can act as an arbitrary phase 



mask, thereby offering a unique platform for simultaneous gain and pulse shaping. This 
shaping is provided by the QPM phase (f)QPM [figure [^f)], which is imparted to the generated 
wave. The idler spectral phase (piiiy) can be approximated as 

oc -0s(z/p -u) + (()QPM{xi{u), 0) - ki{u)L, (3) 

where ki{v) is the idler wavevector and (t)j{v) is the spectral phase of wave j. A more general 
expression for including additional terms to account for the longitudinal variation of 

the QPM grating, is given in Supplementary section 3. 

There is significant freedom in choosing 0), subject only to constraints on the 

ferroelectric domain angles that can be fabricated. Consequently, very large phases can be 
imparted onto the idler. Unlike conventional pulse shapers, this phase is fully continuous, 
and no discrete wrapping between 0 and 27r phase is needed. These properties are illustrated 
in figure [^f), which shows the group delay spectra of the signal and idler for two cases. The 
signal group delay is determined by the dispersion of the material, while the idler group delay 
is determined by equation (|^. The two cases shown correspond to two orientations of a par¬ 
ticular QPM grating, thereby emphasizing how the group delay can be modified substantially, 
over several picoseconds or more, without altering the amplification characteristics. 


4 Experiment 


To demonstrate the technique, we implemented the 2D-QPM FOPA as the final stage of a 


mid-infrared OPCPA system 27,^ (see Methods). A schematic of the FOPA is shown in 


figure l^a). The diffraction gratings have 75 lines/mm, designed for a blaze wavelength of 
4000 nm. The focal length of the 4-f arrangement is / = 200 mm. The elliptical pump beam 
is collinearly overlapped with the spatially-chirped mid-infrared signal by a dichroic mirror 
which transmits the signal and reflects the pump. An identical mirror removes the remaining 
pump after the 2D-QPM crystal (not shown). 
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Figure 3: Experimental setup and results, (a) Schematic of the experimental FOPA setup, 
(b) Input and output spectra; the output spectrum is normalized, while the seed spectrum 
is scaled so that it is visible on the same scale as the pump, (c) Measured and retrieved 
SHG-FROG spectrograms, (d) Reconstructed pulse prohle. (e) Reconstructed spectrum and 
phase. 


The pump laser has a duration of ~ 14 ps, and an average power of 16.5 W at a repetition 
rate of 50 kHz. The beam 1/e^ full-width in the QPM crystal is ^27.4 mm in the horizontal 
and ~ 140 pm in the vertical. The mid-infrared seed has an average power of 6.4 mW before 
the hrst diffraction grating. Its spectral components are spatially chirped according to hgure 
[^b). To improve the temporal overlap with the pump pulses according to equation (|^, we 
use a cylindrical telescope prior to the 4-f arrangement to obtain a relatively large 1/e^ beam 
size of ~ 12 mm (1/e^ full width) in the horizontal direction. 

Measured spectra are shown in hgure |^b). The compressed output power after the second 
diffraction grating was 1.03 W, corresponding to 20.6 pJ pulse energy. Accounting for losses 
of the diffraction grating (?«31 %), the dichroic mirror (~5%), and beam-routing optics, we 
estimate an average power of ~ 1.65 W directly after the antirehection-coated 2D-QPM crys- 
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tal, corresponding to 33 pJ pulse energy. We thus infer a quantum efficiency (ratio between 
output signal photons and input pump photons) of 32%, which is a substantial improve¬ 
ment over our previous OP CPA conhguration based on non-collinear power amplihcation in 


a conventional PPLN crystal 28 


To compress the output pulses, we adjusted the 1550-nm seed pulses in the OPCPA front- 
end. The compressed pulses were measured using second-harmonic generation frequency 
resolved optical gating. The measured and retrieved spectrograms [hgure |^c)] exhibit good 
agreement. The reconstructed pulse prohle is plotted in hgure |^c), indicating compression to 
53 fs (FWHM); this is mainly limited by the available seed bandwidth (43 fs transform limit, 
corresponding to four optical cycles). The reconstructed spectrum and phase are shown in 
figure l^e), in good agreement with the independently-measured spectrum [hgure |^b)]. The 
huctuations on the spectra can be explained by considering the spectral broadening of our 
1550-nm seed in the OPCPA front-end (see Supplementary section 5). 

Our proof of principle experimental results establish the viability of the 2D-QPM FOPA 
technique. The demonstrated conversion efficiency already exceeds the state of the art for 


mid-infrared systems 28 -32 . Nonetheless, the comparatively short duration of the signal 
in the Fourier plane limited the achievable efficiency somewhat. Based on equation (|^, 
using Win{X) = 12 mm (1/e^ full-width), we estimate the signal pulse duration as 6.35 ps 
(full-width at half-maximum), which should be compared to the 14 ps pump. As shown in 
hgure [^d), this mismatch in pulse durations reduces the efficiency compared to the simplihed 
continuous wave case of hgure [^b). Note, however, that since no poling errors were present 
in the fabricated device [hgure [^h)], QPM gratings with signihcantly larger widths up to 
~60 mm are feasible (limited by wafer width). Scaling Ax in equation ([^ accordingly would 
yield a corresponding increase in signal pulse duration. 
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— - - FOPA-1 idler pulse 


Figure 4: Concept for simultaneous amplification and pulse shaping of the input signal 
pulse. Two FOPA crystals are used. The idler generated in the hrst crystal seeds the second 
crystal. In this way, the output of the second crystal inherits the QPM phases of FOPA-1 
and FOPA-2. These QPM phases can be engineered to accomplish the desired pulse shaping 
functionality, for example compensation of higher order spectral phase terms of the signal. 


5 Dual-grating concept for signal shaping 


In this section we propose an extension of the demonstrated approach that enables arbitrary 
shaping of the signal pulse. This is an important point because, unlike the signal, the spatial 
chirp of the generated idler is not linear in wavelength, and hence cannot be fully removed 
by a diffraction grating. The key concept, depicted in figure is to have a second crystal 
seeded by the idler wave from the first crystal. 

We assume an initial signal input to crystal FOPA-1 with center frequency v. In FOPA- 
1, an idler wave (frequency Vp — v) is generated. The spectral phase of this idler is related 
to the QPM phase of FOPA-1 according to equation ([^. Next, the FOPA-1 output signal 
is discarded, while the FOPA-1 idler (frequency Vp — v) seeds FOPA-2. The idler wave of 
FOPA-2 is at the original frequency and its phase is related to the QPM phase of FOPA-2. 
Therefore, the hnal signal output at frequency v of the crystal pair has experienced a phase 
modulation according to the difference of the QPM phases of the two crystals (derived in 
Supplementary section 3). By adjusting the QPM phases by design of the crystals, an arbi¬ 
trary effective phase mask for the signal wave can be achieved, thereby combining amplitude 
and phase shaping into a single crystal-pair. 
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6 Conclusions 


In conclusion, by combining a spatially chirped input wave with a two-dimensionally pat¬ 
terned QPM crystal, we have introduced and demonstrated a new platform for nonlinear 
optics that has unprecedented flexibility and overcomes the limitations and trade-offs inher¬ 
ent in conventional devices. By using QPM crystals with curved domains fabricated with 
high hdelity, we have shown that the technique is scalable in bandwidth, since the QPM 
period trajectory in the crystal can be matched to the input wave, even for extremely broad 
bandwidths where the QPM period changes signihcantly and nonlinearly versus position. 
The approach is applicable to a wide variety of nonlinear-optical devices, including harmonic 
generation and optical parametric amplihcation. 

In contrast to conventional ultrafast processes, the 2D-QPM FOPA consists of a con¬ 
tinuum of narrow-band OPA interactions across the transverse dimension of the crystal, 
with the freedom to adjust the character of these interactions across the spectrum via the 
QPM pattern. For our proof of principle experiment, we combined a transverse variation in 
QPM period with a longitudinal variation to compensate for the pump beam shape, thereby 
addressing two ubiquitous problems in nonlinear optics (phase-matching, and use of non-flat¬ 
top beams). Moreover, a broad class of such longitudinal variations are possible: examples 
include introducing a longitudinal chirp in the QPM period, or introducing a second QPM 
segment for additional functionality such as harmonic generation. Further exploration of 
these capabilities, as well as optimization of the experimental parameters, provides a rich 


framework for future work. Pump beam shaping, as in 14 , could help enable such explo¬ 
ration by reducing constraints due to the pump spatial prohle. 

Beyond amplitude shaping and bandwidth scaling, we showed in section]^ how the FOPA 
can act as an arbitrary phase mask as well, thereby combining the key functionalities of 
broadband amplihcation and pulse shaping in a single crystal-pair. This capability should 
help greatly simplify nonlinear-optical systems, since all of the challenging phase-matching, 
amplihcation, and dispersion management aspects can be lumped into a single, highly engi- 
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neerable component. The slab-like geometry of the FOPA makes it power-scalable, since it 
supports one-dimensional heat flow along the thin axis of the crystal. The advent of large- 


aperture QPM crystals provides great promise for energy scaling as well 33 35 . Moreover, 


QPM media are available in diverse spectral ranges, for example covering from the ultra¬ 


violet via LBGO 36 , to the far-infrared via orientation patterned GaAs and GaP 37 38 


Therefore 2D-QPM devices will be applicable for pulse generation, shaping, and amplification 
across the optical spectrum, from the deep-ultraviolet to far-infrared. We therefore expect 
that frequency-domain processes enabled by such 2D-QPM media will have broad impact 
and appeal for many areas of photonics. 


Methods 

The 2D-QPM grating introduced in section was designed by combining several techniques. 
To obtain the transverse variation of the QPM period (direction x), we calculate the posi¬ 
tion of the spatially chirped spectral components using the grating equation, and apply the 


Sellmeier relation to hnd material phase-mismatch Afco = kp — kg — ki 16 according to hgure 


Bb). For the longitudinal variation (direction z), we first determine the required effective 
length according to figure [^c), and then apply a z-dependent offset in Kg; we construct the 
offset using hyperbolic tangent functions, in analogy to apodization prohles discussed in the 


context of chirped QPM media 20,39,40 . The QPM prohle for a particular x-slice is shown 
in hgure l^a), black curve. 

Rather than simply removing the QPM grating entirely in the switched-off regions, we 
maintain a 50-% QPM duty cycle through the whole crystal, and instead switch off the 
amplihcation process by rapidly changing the QPM period away from phase-matching (i.e. 
a “nonlinear chirp” prohle). Such a 50-% duty cycle helps suppress photorefractive ehects 


41 43 , in addition to the suppression already provided by MgO-doping of the crystal. 


The appropriateness of the chosen QPM designs was tested with simulations of the OPA 
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process including pump depletion. The value of Kg at the end of the grating is not directly 
coupled to the OPA gain provided the phase-mismatch Afc = kp — kg — ki — Kg is made large 
enough. However, nonlinear phase shifts can be introduced by phase-mismatched nonlinear 
processes 


44]: 2D-QPM creates the opportunity to select the phase-mismatch to manipulate 

Aield phase shifts 
. This approach, 


these nonlinear phase shifts. Accordingly, we selected the Kg prohle to yield phase shifts 
comparable in magnitude to those from the intrinsic nonlinearity 
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motivated by our recent work on adiabatic excitation of quadratic solitons 40 , is another 
unique capability of 2D-QPM media which we will explore in more detail in future work. 

With respect to fabrication, our choice of (j)QPMix, 0) in equation Q minimizes the angles 
of the ferroelectric domains. These domain angles satisfy tan(6*(a;, ?/)) = {d(f)QPM/dx)/{d(f)QPM/dz), 
and (1)qpm{x,0) allows the prohle of these domain angles to be manipulated without alter¬ 
ing the longitudinal derivative d(j)QPM/dz which determines the phase-matching properties. 
Lastly, given (1)qpm, the lengths of the pattern on the lithography mask are offset from the 
desired 50-% duty cycle QPM grating to account for domain spreading during poling of the 
ferroelectric domains. 

For our numerical simulations, shown in hgure we developed a unidirectional coupled- 
envelope model designed to model nonlinear mixing between the envelopes due to second- 
and third-order nonlinearities. The model allows in principle for arbitrary mixing between 
the nominal pump, signal, and idler envelopes as well as additional envelopes corresponding 
to sum- and difference-frequencies between the nominal envelopes. The model also supports 
arbitrary QPM media, provided the assumption of unidirectionality and paraxial diffraction 
hold. The spatial chirp of the broadband signal, and the full 2D-QPM phase (j)QPM{x, z), 
are included. The simulation in hgure j^c) uses a “2-l-lD” model (transverse coordinate x, 
propagation coordinate z, and time t). Further discussion is given in Supplementary section 
4. 

Our experiments described in section]^ use a mid-infrared OPCPA front-end containing 


two OPCPA pre-amplihers 27,28 . A schematic of the system is shown in Supplementary 
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figure 1. The system uses two synchronized lasers for pumping and seeding the pre-amplihers. 
The pump laser has a wavelength of 1064 nm and produces 14-ps pulses (FWHM) at a 50-kHz 
repetition rate, with 8 W average power. Approximately 5 W is directed to the pre-amplihers, 
while the remaining power is directed to a home-built Innoslab-type ampliher, the output of 
which is used to pump the frequency domain OPA illustrated in hgurej^a). The seed laser 
is a femtosecond hber laser with subsequent erbium-doped hber amplihers. The laser has a 
wavelength of 1550 nm and produces 70-fs pulses at an 82 MHz repetition rate, with 250-mW 
average power. The seed pulses are hrst spectrally broadened in a dispersion shifted hber 
(DCFS, Thorlabs) before being chirped in time with a silicon prism pair and 4-f pulse shaper 
arrangement. This chirp is transferred to the mid-infrared by the second pre-ampliher, and as 
such we can optimize the compression of the hnal amplihed mid-infrared pulses by adjusting 
the dispersion of the infrared seed pulses. 

The OPCPA pre-amplihers are based on longitudinally chirped quasi-phase-matching 
gratings, implemented in MgOiLiNbOs. We use the shorthand aperiodically poled lithium 
niobate (APPLN) to refer to them in [^. Both crystals have the same QPM design, and 
the pump, signal, and idler beams are all collinearly aligned in the crystals. After the hrst 
APPLN crystal, we discard the idler wave (mid-infrared output), and route the pump and 
signal outputs to the second crystal. After the second APPLN crystal, the pump and signal 
waves are discarded, and we extract the 3400-nm mid-infrared wave to seed the hnal ampliher 
(the FOPA). Further information on the system is given in Supplementary section 1. 

The chirped QPM gratings used for the pre-amplihers can in principle be scaled in band¬ 
width, but careful consideration must be given to several design constraints, described in 


detail in 26 . These constraints, which relate to favoring the desired OPA process over var¬ 
ious unwanted processes, become restrictive when operating in the highly pump depleted 
regime corresponding to adiabatic frequency conversion. Therefore, the combination of lon¬ 
gitudinally chirped QPM devices for convenient and alignment-insensitive pre-amplihcation 
to moderate energy levels, followed by the 2D-QPM FOPA for hnal power amplihcation. 


16 



represents a compelling system arrangement which preserves bandwidth scalability, avoids 
the challenging parasitic processes of highly-saturated chirped QPM devices, and keeps com¬ 
plexity at a minimum since only one FOPA arrangement is required. 
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